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Abstract 

We revisit critically the recent claims, inspired by quantum optics and 
quantum information, that there is entanglement in the biological pigment 
protein complexes, and that it is responsible for the high transport effi- 
ciency. While unexpectedly long coherence times were experimentally dem- 
onstrated, the existence of entanglement is, at the moment, a purely the- 
oretical conjecture; it is this conjecture that we analyze. As demonstra- 
ted by a toy model, a similar transport phenomenology can be obtained 
without generating entanglement. Furthermore, we also argue that even 
if entanglement does exist, it is purely incidental and seems to plays no 
essential role for the transport efficiency. We emphasize that our paper is 
not a proof that entanglement does not exist in light-harvesting complexes 
- this would require a knowledge of the system and its parameters well 
beyond the state of the art. Rather, we present a counter-example to the 
recent claims of entanglement, showing that the arguments, as they stand 
at the moment, are not sufficiently justified and hence cannot be taken as 
proof for the existence of entanglement, let alone of its essential role, in the 
excitation transport. 

1 Introduction 

In recent years, following the development of quantum information, the phe- 
nomenon of quantum entanglement has been identified as being one of the 
most important aspects of quantum mechanics. It was realized that the pres- 
ence of entanglement confers quantum systems significantly enhanced power 
for accomplishing many tasks 12, such as exponentially increased speed-up 
of computation and significantly enhanced communication capacity. As such 
it has been very natural to enquire whether biological systems could have 
evolved to make use of entanglement. Recently the importance of investigat- 
ing this question received a major impetus following seminal experiments that 
indicated the existence of unexpectedly long-time coherent effects in photosyn- 
thesis |2-4|. Since coherence is a pre-requisite for entanglement, its discovery 
raises the possibility that entanglement is also present. Moreover, it is known 
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that the energy transport in light-harvesting complexes is extremely efficient 
- could it be the case that entanglement is responsible for this efficiency? A 
few ground-breaking studies |6j, followed by an increasing body of litera- 
ture JZrlEl raised this question and suggested that this is the case. Here we 
take a critical look at these results. 

The main issue to be discussed here is that of coherence versus entangle- 
ment. They are definitely not one and the same thing. The existence of entan- 
glement is a stronger criterion than the presence of coherence. That is, entan- 
glement requires the presence coherence, but coherence does not imply entan- 
glement in general. While coherent phenomena can also appear in systems of 
classical wave mechanics, entanglement is a genuine quantum phenomenon, 
which is required to violate a Bell inequality, for example. 

It is important to mention from the beginning that while the existence of 
coherence in light-harvesting systems has been experimentally tested (at least 
in laboratory conditions), the existence of entanglement has been not. This is 
not surprising - experimentally proving entanglement is a far more difficult 
task 1131 . As such all the discussions about entanglement in light-harvesting 
systems is, at present, purely theoretical. 

The existence of coherence in these systems seems by now to be well es- 
tablished, and we are not challenging that. It is also quite reasonable to expect 
that coherence plays an important role in the transport problem, distinguishing 
it from classical diffusive processes; we are not challenging this either. What 
concerns us here is the existence of entanglement and its role, if it exists. Specif- 
ically, the main questions we address in our paper are: 

• Are the assumptions that led to the present claims of entanglement in 
photosynthesis justified? 

• Even if these assumptions were justified and entanglement would exist 
along the lines of those models, does the entanglement play any signif- 
icant role in enhancing the efficiency of transport or it is of no conse- 
quence? 

We suggest that, despite the sizable body of literature claiming the contrary, 
the answer to both questions is "no". 

To be clear, we do not prove that there is no entanglement in photosynthe- 
sis. To do that would require knowledge of the system that is well beyond 
what is available at present. All we do here is to present arguments that point 
to potential problems in the present claims that entanglement exists and plays 
a significant role. Our paper should rather be viewed as a counter-example 
and it is meant to sharpen the further investigation of the problem. 

2 Basic Issues 

Light harvesting complexes in plants and photosynthetic bacteria are com- 
prised of protein scaffolds into which pigment molecules are embedded, e.g. 
chlorophyll or bacterio-chlorophyll molecules. The pigment molecules absorb 
light in the visible or infrared part of the spectrum, and the resulting elec- 
tronic excitation (exciton) is transported between the pigment molecules until 
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it reaches a reaction center complex, where its energy is converted into sepa- 
rated charges. 

The FMO protein complex of green sulfur bacteria is a trimeric complex 
that links the chlorosome antenna with the reaction center. Within each of the 
subunits there are seven chlorophyll molecules in close connection with each 
other. Each of these molecules can be considered as a "site" where the excita- 
tion propagating from the chlorosome to the reaction center may be localized. 
It is the entanglement between these sites during the propagation of excita- 
tion that is discussed in 15 HTT11 . (Alternatively so called "mode"-entanglement 
has also been considered in [9|. There, the physical systems between which 
entanglement is investigated are no longer identical to the pigment cofactors 
but effective systems that are defined by the single excitation spectrum of the 
Hamiltonian of the coupled pigments.) 

Regarding the existence and generation of entanglement in light-harvesting 
complexes, we draw some intuition from the following formal analogy. As we 
will detail later, up to certain extent, the FMO complex can be seen analogous 
to a multi-armed interferometer, where each interferometer arm corresponds to 
a site in the FMO complex. The propagation of a single excitation through the 
FMO complex is then analogous to the propagation of a single photon through 
the interferometer. A single photon propagating through an interferometer can 
indeed immediately lead to entanglement between the arms ITHI , hence from 
this point of view it is not very surprising that a single excitation propagating 
through the FMO complex may lead to entanglement between sites. However, 
and this is one of our main concerns, if instead of a single photon we send a 
coherent state through the interferometer, then no arm entanglement will ap- 
pear at all. Indeed, at a beam splitter a coherent state is split into a product of 
coherent states in each of the outgoing arms. As the light passes the various 
beam splitters in the interferometer, in each step we maintain a direct prod- 
uct of coherent states in each of the arms. We stress that this is true even for 
very weak coherent states where the probability of having more than a single 
photon is overwhelmingly small. 

There are two lessons to be learned from this analogy. First, whether or not 
entanglement between interferometer arms exists depends crucially on the ini- 
tial state. Our concern is that in the case of photosynthesis in which the entire 
light-harvesting complex is illuminated by weak classical light and not single 
photons, entanglement may therefore not appear inside the FMO complex. 

The second lesson is that the actual dynamics of light propagation through 
an interferometer, and hence a measure of the transport efficiency, does not re- 
ally depend on the fact that at single photon level entanglement between arms 
is produced. Indeed, the dynamics of an interferometer can very well be de- 
scribed at classical level, i.e. via coherent states, where the question of entan- 
glement does not appear. Hence we conclude that entanglement beyond the 
mere existence of coherences, although it does appear at single photon level, is 
irrelevant for the light transport in interferometers. 

However, excitation transport in the FMO complex is not exactly identi- 
cal to light propagation in an interferometer. It is therefore important to see 
whether or not the problems we mentioned above are still relevant for exci- 
tation propagation through light-harvesting complexes. In this paper we ar- 
gue that, indeed, entanglement plays no essential role in transport through the 
FMO complex. 
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Following from the discussion above, we thus need to revisit the general 
assumptions that have been made in existing theoretical treatments of entan- 
glement in light-harvesting complexes, in particular the following points: 

1. the initial state entering the complex, 

2. the transport between sites and the possible entanglement generated thereby. 

More specifically, we also reconsider the modeling of each of the sites as a 
simplified two-level system. 

3 Initial excitation 

The available literature on entanglement generated during excitation transport 
in pigment protein complexes usually assumes an initial state where (a) only 
one site is excited and (b) this site is excited with exactly one quantum of elec- 
tronic excitation, i.e. a "Fock state", as done in refs. ||5j|6l|8l, for example. Given 
these assumptions, the pigments of the light-harvesting complex have been 
modeled by two-level systems, corresponding to presence or absence of an ex- 
citation at that location. 

The rationale behind the above assumption for the initial state is often the 
following. Under illumination by sun light or, in experiment, by femtosecond 
laser pulses, the photon flux per pigment is weak. An estimation for conditions 
in full sunlight yields a value of about 10 absorbed photons per chlorophyll 
molecule per second in the relevant part of the spectrum Ifl5] . Considering 
that the typical timescales for exciton transport is of the order of picoseconds 
with exciton lifetimes of nanoseconds, this suggests that if light is absorbed, 
then most likely only one excitation is present. 

However, our main point is that the light that excites the light-harvesting 
apparatus is essentially classical, technically a mixture of coherent states. Im- 
portantly even if a coherent state is very weak, and the probability of contain- 
ing more than one excitation is very low, this by no means implies that it can 
be described as a single excitation Fock state. Hence, using the weakness of the 
incident light as an argument for modeling the initial state as a Fock state is a 
fallacy. 

A coherent state differs from a Fock state in two important aspects. First 
of all, a coherent state contains not only one excitation but also terms with 
more excitations. Second, the coherent state also contains a term, the vacuum, 
without excitations. In particular, the existence of the vacuum term has inter- 
esting consequences. It is completely irrelevant as far as the dynamics of exci- 
tation transport is concerned, it just represents the fact that nothing happens, 
but it has crucial implications to entanglement. Since the entanglement that is 
considered in light-harvesting systems stems from superpositions of states in 
which one pigment is excited and all others are in the ground state, a contribu- 
tion with all pigments being in the ground state, corresponding to the vacuum 
in the incoming light, plays an essential role and must not be neglected. 

At this point another possibility arises, namely, that the initial state in the 
FMO complex that is of interest for us, could be a Fock state for other reasons 
than the weakness of the incident light. Indeed, the FMO complex in vivo is 
not directly excited by the incident light, rather it receives its input through the 
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antenna complex. It it thus conceivable that some mechanism in the antenna 
and its connection to the FMO complex may somehow generate a Fock state. In 
quantum mechanical terms this amounts to a state preparation process. What 
would this require? 

Producing a Fock state in the FMO complex, would require two things. On 
one hand, cutting the possibility of more than one excitation. Such mechanisms 
have been discussed 1 16 1, although it is not clear to us how efficient they are in 
limiting the number of higher excitations within the time frame that is relevant 
for transport through light-harvesting complexes. On the other hand, to pro- 
duce a Fock state, one should also be able to eliminate the vacuum component 
of the state. To our knowledge, no such mechanism has been suggested so far. 

This entire discussion above his highly relevant since, as far as entangle- 
ment is concerned, the difference between coherent states, even if they are ex- 
tremely weak, and one-photon Fock states is dramatic. While one-photon Fock 
states may result in entanglement between linearly coupled systems, coherent 
states will not. Furthermore, even in non-linearly coupled systems, weak co- 
herent states can only produce very limit amounts of entanglement, far below 
that produced by one-photon Fock states. 

4 Harmonic oscillator model 

Having discussed the difference between the initial states, we now explicitly 
illustrate the difference that the two states make regarding entanglement. 

In existing studies of entanglement in light-harvesting complexes, the man- 
ifold of electronic states of each of the relevant pigments, e.g. chlorophyll mol- 
ecules, is usually restricted to only a few levels, mostly only the electronic 
ground state and the first excited state. The rich energetic landscape of each 
chlorophyll molecule is thereby conceptually replaced by a two-level atom. 

In order to account for small contributions of higher excited states, we must 
consider a different model for the electronic level structure of each pigment 
molecule, that allows for more excitations to be present at a single site. As the 
simplest possible scenario, consider the analogy of a light-harvesting complex 
to an interferometer. In this analogy, the electronic level structure of each site 
is modeled by a harmonic oscillator and thus it is formally identical to a light 
mode in an interferometer. 

For capturing the principles of entanglement generation during the evolu- 
tion of the excitations in a network of coupled chlorophyll pigments, it is suffi- 
cient to first study the interaction and state evolution of the simplest network 
of only two sites, a dimer. Hence, we first consider how entanglement is gener- 
ated between two coupled harmonic oscillators. (Incidentally, the present ap- 
proach to model the pigment molecules by harmonic oscillators has also been 
employed in [17|, however, by means of classical harmonic oscillators and with 
a different aim, namely to illustrate that the phenomenology of quantum co- 
herent exciton transport can also be obtained with a classical coherent model. 
In the present work, however, we employ a full quantum description, and our 
focus lies on entanglement.) 

We assume for the interaction Hamiltonian a standard form where excita- 
tions are exchanged between the modes, the rotating -wave approximation has 
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already been applied, and the systems are taken to be resonant: 

Hin t =^(fl + b + ^ + ). (1) 

The coupling strength is denoted by g, and the creation and annihilation oper- 
ators a, b, and a + , for the harmonic oscillators describing molecules A and 
B, respectively realize the exchange of an excitation between the molecules. 

To exemplify our argument with coherent states, we investigate the initial 
state of a coherent state for molecule A, which expanded into the Fock basis of 
states with n excitations reads 

hfc(0)) = |*)=e-N /2 E^I">' ( 2 ) 

and whose mean number of excitations is given by \oc\ 2 . For molecule B we 
take the vacuum (ground) state \Tps(0)) = |0). The requirement that the light 
intensity is low and hence an excitation occurs only with little probability, and 
that higher excited states should virtually not occur, thus formally amounts to 
\a\ < 1 for the initial state. This initial state is in accordance with the cen- 
tral assumption that if an excitation occurs, then most probably there is only a 
single excitation, because for n > 1 

|(1|^(0)>| 2 = e-H 2 | a |2 » e-^ 2 \K\ 2 "/n\ = \(n\ip A (0))\ 2 . (3) 

In the interferometer analogy this initial state amounts to a coherent state and 
a vacuum state for the two incident light modes, respectively. 

Under the given Hamiltonian, the state of the two molecules evolves in the 
interaction picture according to 

\xp AB (t)} = U(t)\ip A (0)}\MO))=e- lsHabf+a * h) l»)|0). (4) 

The time evolution of this system of two coupled harmonic oscillators is an 
elementary problem. A standard textbook calculation yields ITT8I 

\rp(t)) = U(t)D A (a)U\t)U(t)\0)\0) (5) 

= e ^U^a*UaU f |Q> |0>. (6) 

Here, the displacement operator D^(a) = exp(aa + — m*a) is used to construct 
coherent states by simply moving the ground state (vacuum) away from the 
phase space origin, \a) = Da(#)|0). The unitarity of U(t) allows to sandwich 
all operators in the exponential by inserting identity operations, the time argu- 
ment has been omitted, and the ground state is invariant under U. Using the 
relation, 

e x Y e- x = Y + [X, Y] + 1 [X, [X, Y}] + 1 [X, [X, [X, Y]]] + ■ ■ ■ , (7) 

the sandwiched creation and annihilation operators can be evaluated: 

= cos(#f > + + i sin(gf )k + (8) 
UaU f = cos{gt)a - i sisx(gt)b. (9) 
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Since the operators for molecule A and B commute, the exponential can be 
split in two parts that each give a displacement operator for molecule A and B, 
respectively, 

\ip(t)) =D A [acos(gt)] D B [iasin(gt)] |0>|0) 

= \occos(gt)) |zasin(gf)). (10) 

In the final state after some interaction time t, each of the two molecules is in 
a coherent state with parameters acos(gf) and iasin(gt), respectively. With 
respect to entanglement, let us point out that for all interaction times and ir- 
respective of a, the two molecules are always in a product state, i.e. there is 
strictly no entanglement between the molecules. 

Let us spell out the formal analogy between the dynamics in our simplified 
model of the two interacting chlorophyll molecules, the electronic structure of 
each of which is modeled by a harmonic oscillator, and two light modes inter- 
acting at a beam splitter. At a beam splitter, in Schrodinger picture, the state of 
the two incident and outgoing light modes \<pi n ) and \(p ou t), respectively, is re- 
lated by a unitary operation 1 19 1, which contains the reflectivity /transmissivity 
of the beam splitter as a parameter 6: 

Ifout) = liWIfcn) = e-^ 8 ^) (11) 

The unitary operation of the beam splitter is formally identical to that of the 
two interacting molecules in Q. Therefore, the state of the molecules before 
and after they have interacted for some time t can be associated to the state of 
the two light modes before and after they have passed the beam splitter with 
reflectivity/ transmissivity parameter 8/2 = gt. A 50/50 beam splitter with 
9 = 7i / 2 thus gives states in the output ports that correspond to the states of 
the molecules after time gt = tt/4: For an initial coherent state, one obtains a 
product of coherent states as in (TO) . In contrast, for a single photon entering 
the beam splitter in one input port, \(pi n ) = |10), the state emerging from the 
beam splitter is, |^ ut) = ( 1 10} + z 1 1 } ) /y/2, a maximally entangled state of the 
two interferometer arms. 

Let us translate the observation of a coherent state versus a single photon 
Fock state from the interferometer analogy back to the model system of inter- 
acting pigment molecules. Given that the initial state is a coherent state, the act 
of limiting the focus of the treatment to the single excitation manifold creates 
the illusion of entanglement. Formally, one projects the total state | (10) to the 
subspace with exactly a single excitation and obtains for the initial state 

i\|«>|0)«|10), (12) 

where 

Pi = |10)(10| + |01)(01| (13) 

is the projection onto the single excitation subspace. That is, there is one exci- 
tation on molecule A and none on B. For the time-evolved state this projection 
yields 

P\\$ab (0> = Pi|«cos(gf)) | fa sm(gt)} ck cos (gt) 1 10) + i sin(^)|01). (14) 
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In particular for gt = n/A the state appears to be maximally entangled in anal- 
ogy to a single photon traversing a 50 /50 beam splitter. 



It is interesting to ask why does considering all possible excitations remove 
the entanglement? How can the higher excitations sector, which has a very 
small contribution, remove the entanglement? This seems to be quite paradox- 
ical. The answer is that if the coherent state is so weak as to have only a very 
small amount of higher excitations, it necessarily also has a very large vacuum 
component. Just considering the vacuum term in addition to the single exci- 
tation sector reduces the amount of entanglement considerably. Any residual 
entanglement is eliminated by the higher excitation sector. 

To see the effects from above, it is convenient to use a measure for the 
amount of entanglement as given by the concurrence |20], which for a pure 
state of two systems is defined in terms of the purity of one (any) of the sub- 
systems: 



C(Hab}) = ^2(1-Trp2) W i t h PA =Tr B \y AB ){ip AB \. (15) 

The concurrence of the (renormalized) state after projecting out the single ex- 
citation sector gives 

c(l£^)=\M2gt)\, (16) 



AB/ 



which is maximal for gt = tz/4. 

Alternatively, if the attention is not strictly limited to the single excitation 
manifold but the vacuum (ground state) term is considered in addition, the 
relevant projection to apply to the state is 

P 0< i = |00)(00|+Pi. (17) 

For the initial state, this projection yields a superposition of the ground and 
excited state of the first site, i.e. 

%|tMB(0)) = Po,lk>|0) « (|0)+«|1))|0), (18) 

and, for the time-evolved state, a superposition of the ground state and the 
evolved state in the single excitation manifold: 



PoaWM*)) k l°°> + *[a»(gt)\l0) +isin($t)|01)J. (19) 
The entanglement in terms of concurrence of the normalized state is therefore 



\Po,i\iPab) 

We find that the maximal amount of entanglement at gt = 7i/4 is limited by 
the dominant ground state contribution for small \ct\, because the ground state 
is not entangled. 

The entanglement that appears after projecting the state onto its zero and 
single-excitation subspace with Pq \ is simply rescaled with respect to the en- 
tanglement found when projecting only to the single excitation manifold. The 
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scaling factor |a| 2 /(l + |a| 2 ) amounts to the single excitation fraction of the 
projected state, i.e. the probability of measuring exactly a single excitation in 
the state once it has been projected. As long as the assumption holds that the 
probability of more than one excitation is small, we require that also | a\ is very 
small, and therefore the entanglement reduction due to ground state contribu- 
tion is larger. The higher excitation terms in the coherent state may be small in 
absolute magnitude, but they only have to eliminate this weak residual entan- 
glement. 

Another aspect that needs to be taken into account when advancing to more 
realistic models for the study of quantum entanglement in light-harvesting 
complexes is that superpositions of electronic states on the same molecule may 
suffer from decoherence, primarily due to the coupling of the electronic structure 
to nuclear degrees of freedom. This process ultimately turns a coherent state 

p A (0) = |^(0)>(^(0)| = |*X*| / (21) 

when written as a density operator, into the incoherent mixture 

n 00 U|2n 

pA (t^oo)=e-W EMH"X"I (22) 

n=0 

of electronic states of the individual molecule, meaning that after complete 
decoherence all coherences between states of different excitation number will 
have decayed. At room temperature molecules have been demonstrated to re- 
main in a coherent superposition of their electronic ground and first excited 
state for a 50 fs-timescale |21j, which might, however, differ for pigment mole- 
cules embedded in a protein matrix. 

In the worst case, a complete decoherence of an initially coherent state and 
subsequent projection to ground and single excited manifold yields the inco- 
herent mixture of the ground state and the coherently propagated single exci- 
tation: 



|00)(00| (23) 



1 + M 2 

+ |a| 2 (cos(gf)|K)) +isin(gf)|01>) (cos(gf)(10| - zsin(gf)(01| 



The concurrence for this mixed density matrix can be evaluated with Wootter's 
formula l20l , and gives 



c(p^r h \t))- l j ^\M2gt)\, (24) 
that is, the same result as for the projected coherent state l |2u] . 



5 N-level system model 

After we have formally established that the molecules are not entangled if the 
initial state is a coherent state, and that a projection of the state to the single 
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excitation manifold creates the illusion that entanglement were present, let us 
now investigate how the number of considered levels in addition to the ground 
state and the single excited state affects the apparent amount of entanglement. 
In other words, we generalize from the projections Pi and Po,i of the previous 
section to projections that include higher number of excitations. 

An alternative way to look at this question is the following consideration. 
The electronic eigenstates of molecules do not form an infinite uniform level 
structure as a harmonic oscillator. It is reasonable to assume that only finitely 
many excitations can be supported by each pigment molecule, and any num- 
ber of excitations beyond a certain threshold would cause ionization processes 
that take this fraction of the state out of the considered events. For example, 
instead of a coherent state at one of the molecules, which involves contribu- 
tions of arbitrarily many excitations, the molecule may at most support two 
excitations. The initial state would thus be the first three terms of the coherent 
state expansion until n = 2: 

Itf } (0)> = 1 + H ^ k|4/2 , (|0> + «H> + ^|2>) , (25) 

which corresponds to the projection Pqx2 applied to the coherent state \ot), and 
the result renormalized. For these kinds of initial states, which constitute a 
modification of coherent states towards a more realistic initial state for pig- 
ment molecules, we now study how the number of additionally considered 
contributions of higher lying excited states affects the amount of entanglement 
that is generated during the transfer of excitations between two pigment mole- 
cules. Clearly, for the highest considered excited state being N = 1, we recover 
the result | [20) , whereas in the limit N — >■ oo we recover the case of coherent 
states that never generate entanglement. 

In general, we model each molecule by system of N levels, i.e. the ground 
state and N — 1 excited states. In parallel with the harmonic oscillator, we 
choose the initial state to be the ground state for molecule B, and a state for 
molecule A, 

|^ N) (0)> = M = -= £ -=|n), (26) 

which is the projection of a coherent state to the lowest lying N levels. The 
squared norm A/^n = Y^=o\ 0i \ 2 / nl - approaches the value el "1 2 found for co- 
herent states in the limit N — » oo. 

As another exemplary case, let us give the analytic expressions for N = 
3 for the initial state ( (25} as it evolves according to the interaction Hamilto- 
nian Q. The time-evolved state necessarily only contains terms with at most 
two excitations: 

^(0) = v / 1 + k |2 + k |4 7i(l 00 ) +ffi ( COS( ^ )|10) ~ fSin( ^ )|0i: 
2 \ 

+ -= (cos 2 (g0|20) -sin 2 (g0|02) -iV2cos(gt)sm(gt)\ll)) J. (27) 

For \a\ <C 1, the doubly excited states give only a perturbative correction to the 
expression of the quantum state to the case N = 2 in < [T9) . The expression for 
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Figure 1: Scaling of maximal value of apparent entanglement as measured by 

concurrence, C^aLr ror different a as a function of the number of considered 
levels N. 



entanglement, however, does not only change by a perturbative correction, but 
it changes considerably. The concurrence of the state is given by 



r /„(3 K \ k| 3 |sin(2^)|V8+2l«l 2 ( 1 3 + 3cos(4^)) 

C U^)J - 4(1 + W 2 + W 4 /2) (28) 

with its maximum at gt = tc/AloI 

(3) _ MVsTW (29) 

max " 4(l + |«|2 + | a |4/ 2 )- 

The expression of the concurrence for N = 3 is similar in structure to that of 
the case N = 2 in | [20) , with a |sin(2g£)| modulation in time, and the squared 
norm of the projected coherent state in the denominator. The time-independent 
prefactor, however, can no longer be interpreted anymore as the probability of 



(2) 

having an excitation in the system, as done for C m i x in 1 20 1. Instead, we find a 
scaling with the third power of \oc\. 

Expressions for larger N can be straightforwardly obtained, but are omitted 
here. We have analytically examined expressions for the concurrence up to 
N = 7 and generally find a global |sin(2gf)| modulation in time in accordance 
with the intuition that maximal entanglement for these initial states is reached 
for the 50/50 beam splitter configuration, that is, after half the time that an 
excitation needs to fully move from one pigment to the other. In particular, the 
expressions for maximal entanglement are of the generic form 

c£2 = J^/n(V| 2 ), (30) 



N 



that is, the apparent amount of entanglement decreases exponentially with the 
number of considered levels per molecule for \u\ < 1 . The factor fa for higher 
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0.7071 


0.3819 


0.1768 


0.0734 


0.0280 



Table 1: Constant leading order terms of fa (\ot\ 2 ) for |a| <C 1. 



N is of a similar square root form as found for N = 3. For |a| < 1, the leading 
term of fa (|a| 2 ) is constant but also decreases with N, as listed in table 1 
Figure [T] shows how the maximal apparent amount of entanglement decreases 
with N for various choices of a. 

Obviously, the interpretation that given an initial coherent state for \oc\ -C 1 
one arrives at the single photon level is misleading because the entanglement 
properties of a single photon Fock state are not obtained in this limit. Instead, 
for smaller values of \oc\, the apparent entanglement will vanish more rapidly 
when increasingly many levels of the system are taken into account although 
they contribute to an ever smaller extent (see figure 111. 

We have also investigated an alternative way of modeling the initial state 
of an N-level molecule by means of atomic coherent states (also called spin co- 
herent states) ||18l l22l , where instead of the raising and lowering operators of 
the harmonic oscillator, those of the angular momentum algebra are employed. 
The N levels of each molecule are thereby modeled by the N = 2s + 1 levels 
of an effective spin-s particle. Although the computations are more intricate 
due to different commutation relations, we arrive at qualitatively similar re- 
sults as presented for the projected harmonic oscillator. We therefore conclude 
that our argument concerning the entanglement content of system that starts 
in a coherent state is robust with respect to the specific framework applied. 

6 Transport efficiency versus entanglement 

The excitation transfer efficiency in chromophore complexes captures how well 
an excitation that starts somewhere localized in the complex traverses the net- 
work of coupled chromophores to a different location, where it is assumed to 
leave the complex, e.g. to the reaction center. A question of current interest 
that has been addressed |8 11 J is, whether or not entanglement (rather than 
mere coherence) impacts the transport efficiency. 

Let us first recollect a few essential facts about the recent treatments of 
excitation propagation and evaluation of transport efficiency in conjunction 
with the study of entanglement in light-harvesting complexes. In a complex of 
coupled chromophores with pairwise interaction Hamiltonians of the kind as 
in fill, the number of excitation quanta are conserved. Therefore, subspaces of 
a fixed number of excitation quanta, e.g. the single-excitation subspace, evolve 
independently from each other. Cross-contribution of subspaces with a fixed 
excitation number occur only due to interaction of the chromophores with 
other degrees of freedom, i.e. an environment, and for the excitation energies 
considered here (~1 eV) mostly downward to lower numbers of excitation, due 
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to relaxation for example. Given that the assumption holds that the light inten- 
sity is weak and therefore the presence of higher number of excitations hap- 
pens only with minuscule probability the single-excitation subspace is essen- 
tially only subject to the unitary dynamics according to the system Hamilto- 
nian, excitation-number conserving environment influences such as decoher- 
ence, and (excitation-number non-conserving) relaxation to the ground state. 
Therefore, the single excitation-subspace evolves largely independent for the 
other subspaces even in the presence of decay mechanisms, and in particular 
it does not significantly gain excitations from higher lying states during the 
transfer through the complex. Excitation transfer can only occur via the ex- 
cited states, since only then an excitation (at least one) is present. The transfer 
efficiency is usually evaluated by an observable, which is therefore defined 
only in the excited state manifold, and for the scenarios considered here, gains 
its dominant contribution from the transport dynamics of the single-excitation 
manifold. Contributions from higher excited states constitute only a pertur- 
bative corrections to the quantum state and therefore also to the observable 
that quantifies the transport properties. Common examples for measures of 
the transport efficiency are obtained by integrating the population of a certain 
exit site, e.g. decay to the reaction center, as done in |5|, or by the highest pop- 
ulation of an exit-site during a certain time-window as in [11 J, for example. 

Formally, for an observable T that measures the excitation transport effi- 
ciency in the described way, and for initial states with only a small or even 
vanishing fraction of higher excited states, one has 

(T) » (PiTPi), (31) 

that is, one can restrict the evaluation of T with good agreement to the predom- 
inant contribution from the single excitation manifold, because higher excited 
states yield only a perturbative correction to this result. Since the Hamiltonian 
conserves the number of excitation quanta, it is a valid approach to restrict the 
propagation of the entire excitation dynamics to the single-excitation subspace 
right from the beginning, as used in the last step of the following transforma- 
tion: 

(T)« (ip(t)\PiTPim)) 02) 
= (^(0)|U + (f)PiTPiU(0|!/'(0)) (33) 
= (ipWlPirfWTUWPiMO)). (34) 

Because the projection to the single-excitation manifold is effectively carried 
out when measuring the transport efficiency, the projection may as well be 
exchanged with the dynamics to the beginning of the process. 

Even with respect to open system dynamics, which is captured by the dy- 
namical map A(f), that is, the map that contains the formal solution to the 
Liouville equation p(t) = Cp{i), one can approximate the transport efficiency 
by only considering the single excitation manifold and the ground state: 

(T) = Tr [Tp(t)\ = Tr [TA(t)p(0)] » Tr [TA(f)%p(0)P ,i] • (35) 

Regarding the evaluation of a measure of transport efficiency it is thus, for the 
given assumption of weak light intensity, a natural and justified assumption to 
restrict the investigation to the subspace with only a single excitation. 
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In contrast to observables like measures of transport efficiencies, entangle- 
ment must be a non-linear property of quantum states, and it can thus generally 
be not equivalent to coherences, which can be extracted by a linear operator. 
Because of these intrinsic properties of entanglement measures, it is not pos- 
sible to exchange in a similar way the dynamics in the full space with a part 
that is projected to the single-excitation subspace. In fact, it is the projection 
to a subspace of fixed excitation number, which is a global operation, that in- 
troduces the observed entanglement into the system in the first place. In the 
present case of two molecules, the projection Pi acts like a Bell-state measure- 
ment. 

From the presented case we can now provide an insight about whether or 
not the entanglement that may be observed in manifolds of fixed excitation 
number is of relevance to the state evolution or the transport properties of 
the system. The transport efficiency is robust under changes of the underly- 
ing models and assumptions about the initial state regarding the presence of 
small perturbative corrections of higher lying excited states, whereas entangle- 
ment is not. Therefore, in a pigment protein complex where an initial excitation 
merely evolves according to the system Hamiltonian and under a coupling to 
a bath, which introduces decoherence or relaxation, entanglement cannot be 
a quantifier of transfer efficiency. It cannot be the cause of a large transport 
efficiency, nor enhance the transport efficiency. The propagation of an excita- 
tion in a pigment protein complex is different from the case of, say, quantum 
information communication tasks like quantum teleportation, where entangle- 
ment can be identified as the key resource and unentangled (separable) states 
cannot be used or, as another example, the fact that in interacting quantum 
systems entanglement is required to reach the ground state. 

To conclude, in the first place, it is not at all clear to us that entanglement 
exists in the FMO complex; in fact, it seems to us that it does not. Entangle- 
ment, as it is at present postulated, relies on the existence of a single excitation 
in the FMO complex (technically a single excitation Fock state). As we argued, 
this cannot be obtained by simply having very weak light impinging on the 
light harvesting complex, as wrongly assumed in most literature on the sub- 
ject. Consequently, if not simply due to low intensity light, the only other way 
in which the postulated entanglement may still exist is if there would be a 
dynamical, active, mechanism of state preparation. We argue that this is also 
extremely unlikely. Indeed, one would need a non-trivial process based on 
measuring the number of excitations, and opening the entrance of the FMO 
complex depending on the presence of a single excitation. Crucially, this pro- 
cess should also eliminate the large vacuum component. Finally, and more 
significantly, even if entanglement exists, its role for the transport efficiency 
seems to us to be irrelevant, whereas the role of coherence may be important. 
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